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FOREWORD 


This report constitutes a part of the work completed during the year 
1975-1976 on a research project entitled "Radiation-Induced Precursor 
Flow Field Ahead of an Entry Body in the Outer Planetary Atmospheres." 

The work was supported by the NASA-Langley Research Center through 
Contract NAS1-11707-92. The contract was monitored by Mr. Randolph A. 
Graves, Jr. of SSD-Aero thermodynamics Branch. 


ii 


TABLE OF CONTENTS 




Page 


FOREWORD ii 

SUMMARY 1 

LIST OF SYMBOLS 2 

1. INTRODUCTION 5 

2. BASIC FORMULATION 8 

3. SOLUTIONS FOR SPECIAL CASES 15 

3.1 Radiation from a Plane Source 15 

3.2 Radiation from Spherical and Cylindrical Point Sources . 17 

3.3 Solutions for the Transition Range 21 

4. PERTURBATION EQUATIONS FOR PHOTOABSORPTION MODEL .... 22 

4.1 Photoabsorption Model 22 

4.2 Precursor Equations for the Photoabsorption Model . . 24 

5. DATA SOURCE AND SOLUTION PROCEDURE 28 

6. PERTURBATION RESULTS 38 

•7. ALTERNATE APPROACH: THIN-LAYER APPROXIMATION 61 

8. CONCLUSIONS 70 

REFERENCES 72 

APPENDIX Al: EXPLANATION OF SYMBOLS USED IN THE COMPUTER 

PROGRAM "PERC" Al 

APPENDIX A2: EXPLANATION OF SYMBOLS USED IN THE COMPUTER 

PROGRAM '•THIN" A3 

APPENDIX Bl: PROGRAM "PERC" - TO COMPUTE PRECURSOR EFFECTS 

USING PERTURBATION METHOD Bl 

APPENDIX B2: PROGRAM "THIN" - TO COMPUTE PRECURSOR EFFECTS 

USING THIN-LAYER APPROXIMATION B5 


iii 



LIST OF FIGURES 


Figure Page 


2.1 Physical model and coordinate system for a Jovian 

entry body 9 

3.1 Approximation of radiating flow by point sources ... 18 

4.1 Absorption cross section of in ultraviolet region . 23 

5.1 Atmospheric conditions for Jupiter's entry 29 

5.2 Inverse mean free path and (2R C ) as a function of 

altitude for the three regions of photoabsorption 

model 31 

5.3 A(v) as a function of free-stream velocity for spectral 

' ranges I, II, and III (figs, a-c) 33 

6.1 Velocity perturbation as a function of distance from 

the shock at different altitudes and a constant free- 
stream velocity (figs, a and b) 39 

6.2 Pressure perturbation as a function of distance from 

the shock at different altitudes and a constant free- 
stream velocity (figs, a and b) 41 

6.3 Mass fraction of H as a function of distance from 

the shock at different altitudes and a constant free- 
stream velocity (figs, a and b) 43 

6.4 Mass fraction of H 2 + as a function of distance from 

the shock at different altitudes and a constant free- 
stream velocity (figs, a and b) 45 

1 

6.5 Temperature perturbation as a function of distance from 

the shock at different altitudes and a constant free- 
stream velocity (figs, a and b) 47 

6.6 Specific total enthalpy perturbation as a function of 

distance from the shock at different altitudes and a 
constant free-stream velocity (figs, a and b) 49 

6.7 Velocity perturbation (just ahead of the shock) as a 

function of free-stream velocity for constant 

altitudes 52 

6.8 Pressure perturbation (just ahead of the shock) as a 

function of free-stream velocity for constant 

altitudes 53 

(cont'd.) 

iv 



LIST OF FIGURES - CONCLUDED 


6.9 Mass fraction of H (just ahead of the shock) as a 
function of free-stream velocity for constant 
altitudes 


6.10 Mass fraction of H2+ (just ahead of the shock) as a 
function of free-stream velocity for constant 
altitudes 


6.11 Temperature perturbation (just ahead of the shock) as a 
function of free-stream velocity for constant 
altitudes 


6.12 Specific total enthalpy perturbation (just ahead of the 
shock) as a function of free-stream velocity for 
constant altitudes 


6.13 Absolute enthalpy (just ahead of the shock) as a 
function of free-stream velocity for constant 
altitudes . 


6.14 Increase in total enthalpy (just ahead of the shock) as 
a function of free-stream velocity for constant 
altitudes 


7.1 Curvilinear orthogonal coordinate system for thin-layer 
approximation 


7.2 Comparison of results for velocity variation in the 

precursor zone 67 


7.3 Comparison of results for pressure variation in the 
precursor zone 


7.4 Comparison of results for temperature variation in the 

precursor zone 69 


LIST OF TABLES 


Free-stream and shock conditions for Jovian entry 


Temperature coefficients for thermodynamic functions 
for hydrogen species 


v 


1 ■l rg r '»-- 


RADIATION INDUCED PRECURSOR FLOW FIELD 
AHEAD OF A JOVIAN ENTRY BODY 

by 

S. N. Tiwarf* and K. Y. Szema** 

SUMMARY 

The change In flow properties ahead of the bow shock of a Jovian entry 
body, resulting from absorption of radiation from the shock layer, is inves- 
tigated. Ultraviolet radiation is absorbed by the free stream gases, causing 
dissociation, ionization, and an increase in enthalpy of flow ahead of the 
shock wave. As a result of increased fluid enthalpy, the entire flow field 
in the precursor region is perturbed. The variation in flow properties is 
determined by employing the small perturbation technique of classical aero- 
dynamics as well as the thin layer approximation for the preheating zone. 

By employing physically realistic models for radiative transfer, solutions 
are obtained for velocity, pressure, density, temperature, and enthalpy 
variations. The results indicate that the precursor flow effects, in general, 
are greater at higher altitudes. Just ahead of the shock, however, the 
effects are larger at lower altitudes. Pre-heating of the gas significantly 
Increases the static pressure and temperature ahead of the shock for velo- 
cities exceeding 36 km/sec. The agreement between the small perturbation 
and thin layer approximation results are found to be excellent. 


^Professor, Old Dominion University 
**Research Assistant, Old Dominion University 
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pressure perturbation 

free stream pressure, dyne/ cm 2 

divergence of net radiant heat flux, erg/cm 3 -sec 

specific radiative flux density at shock wave, erg/cm 2 

universal gas constant * 8.3143 * 10 7 erg/mole-°K 

radius of the radiating gas cap, cm 

radius of the bow shock wave, cm 

cylindrical radial coordinate, cm 

spherical radial coordinate, cm 

temperature, °K 

temperature at the shock front, °K 
free-stream temperature, °K 
first-order temperature perturbation, °K 
second-order temperature perturbation, °K 
velocity vector, cm/sec 

first-order velocity perturbation, cm/sec 

second-order velocity perturbation, cm/sec 

velocity component normal to the shock surface, cm/sec 

free-stream velocity, cm/sec 

velocity perturbation component at the r-direction in cylindrical 
coordinate, cm/sec 

velocity perturbation component at the x-direction, cm/sec 

velocity perturbation component at the y-direction, cm/sec 

velocity perturbation component at the z-direction, cm/sec 

molecular weight, gm/mole 
photodissociation yield 

photoionization yield 
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Z altitude of entry, km (Table 5.1) 

z longitudinal coordinate, cm 

$ quantity defined in Eq. (4.5) 

r quantity defined in Eq. (2.23a) 

Y specific heat ratio 

e effective emissivity 

C optical depth defined in Eq. (3.3) 

ri quantity defined in Eq. (3.14) 

inverse of photon mean free path, 1/cm 
A quantity defined in Eq. (3.14) 

U quantity defined in Eq. (3.14) 

P density, g/cm 3 

first-order density perturbation 

o Stefan-Boltzmann constant ■ 5.6697 * 10~ 5 erg/cm 2 -sec- 8 K 4 

0p(v) photodissociation absorption cross section, cm 2 

Oj(v) photoionization absorption cross section, cm 2 

4 potential function defined in Eq. (2.27) 

<p potential function defined in Eq. (2.13) 

ip quantity defined in Eq. (4.5) 
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1 . INTRODUCTION 


The word "precursor" gets its name from a Latin word "praecursor" 

(prae * before + currere ■ run) which means "forerunner." j n present 
context, precursor region flow (or flow in the precursor zone) means flow 
field ahead of a shock layer which is influenced by the conditions of the 
shock layer. It is well understood now that at high speed entry conditions 
(entry speeds in excess of parabolic speed) , radiation plays a very important 
role in the analyses of flow phenomena around the body and that the radiative 
energy transferred to the body usually overtakes the aerodynamic heat transfer 
[1-10].* Radiative energy transfer from the shock layer of a blunt body into 
the free stream reduces the total enthalpy of the shock layer while increasing 
the enthalpy of the free stream gases. Because of this increase in enthalpy 
the entire flow field ahead of the shock layer and around the body is influ- 
enced significantly. The precursor flow region is considered to be the region 
ahead of a shock wave in which the flow field parameters have been changed from 
free stream conditions due to absorption of radiation from the incandescent 
shock layer. Most of the radiative energy transferred from the shock layer 
into the cold region ahead of the shock is lost to infinity unless it is equal 
to or greater than the energy required for dissociation of the cold gas. When 
the photon energy is greater than the dissociation energy, it is strongly ab- 
sorbed by the cold gas in the ultraviolet continuum range. The absorbed energy 
dissociates and ionizes the gas and this results in change of flow properties 
in the precursor region. In particular, the temperature and pressure of the 
gas is Increased while velocity is decreased. The change in flow properties 
of the precursor region, in turn, influences the flow characteristics within 
the shock layer itself. The problem, therefore, becomes a coupled one and 
iterative methods are required for its solution. 


*Numbers in brackets Indicate references. 
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Only a limited number of analyses on radiation-induced precursor flow are 
available in the literature. Works available until 1968 are discussed, in 
detail, by Smith [11,12], By employing the linearized theory of aerodynamics. 
Smith [11,12] investigated the flow in the precursor region of a reentry body 
in the earth's atmosphere. The cases of plane, spherical, and cylindrical 
point sources were considered and solutions were obtained for a range of alti- 
tudes and free stream conditions. It was found that for velocities exceeding 
18 km/ sec, precursor flow effects are greatest at altitudes between 30 and 46 km. 
It was further concluded that preheating of air may cause an order of magnitude 
Increase in the static pressure and temperature ahead of the shock wave for 
velocities exceeding 15 km/sec. Lasher and Wilson [13,14] investigated the 
level of precursor absorption and its resultant effect on surface radiation 
heating for earth's entry conditions. They concluded that, for velocities less 
than 18 km/sec, precursor heating effects are relatively unimportant in deter- 
mining the radiative flux reaching the surface. At velocities greater than 
18 km/sec, the amount of energy loss from the shock layer and resultant pre- 
cursor heating correction was found to be significantly large. Liu [15,16] also 
investigated the influence of upstream absorption by cold air on the stagnation 
region shock layer radiation. The thin layer approximation was applied to both 
the shock layer and the preheating zone (the precursor region). The problem 
was formulated for the inviscid flow over smooth blunt bodies but the detailed 
calculations were carried out only for the stagnation region. The general 
results were compared with results of two approximate formulations. The first 
approximate formulation neglects the upstream influence and the second one 
essentially uses the iterative procedure described by Lasher and Wilson [13,14]. 
The results are compared for different values of a radiation/convection para- 
meter. A few other works, related to the effects of upstream absorption by air 
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on Che shock layer radiation, are discussed by Liu [15,16]. Some works on 
precursor ionization for air as well as hydrogen-h( .ium atmosphere are pre- 
sented in [17-21]. 

The purpose of this study is to Investigate the changes in flow properties 
in the preheating zone of a Jovian entry body resulting from absorption of the 
radiation from the shock layer. As a first approach, the perturbation tech- 
nique adapted by Smith [11,12] for the earth's atmosphere is used here for the 
hydrogen-helium atmosphere. By Introducing appropriate thermodynamic and 
spectral information on hydrogen-helium atmosphere, proper modifications are 
made in the governing equations and results are obtained for Jupiter's entry 
conditions. Basic formulation of the problem is presented in Sec. 2 and solu- 
tions for special cases are obtained in Sec. 3. Perturbation equations are 
specialized for the photoabsorption model in Sec. 4, sources of data and solu- 
tion procedures are discussed in Sec. 5, and results of flow perturbations are 
presented in Sec. 6. Finally, in Sec. 7, an alternate approach of thin pre- 
cursor region approximation is adopted for analysis of the entire problem and 
the results are compared with the results of small perturbation theory. 
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2 . BASIC FORMULATION 


The physical model and coordinate system for a Jovian entry body is shown 
in Fig. 2.1. The flow field ahead of the body can be divided primarily into 
two regions, the precursor region and the shock layer. In this study, attention 
is directed to the precursor region where flow is assumed to be steady and invis- 
cid. The flow properties are considered to be uniform at large distances from 
the body. For this region, conservation equations can be written as [22-24] 


Mass Continuity: V • (pV) =0 (2.1) 

Momentum: p(V • VV) = -Vp (2.2) 

Energy: p(V • VH^) - Q R (2.3) 

Species Continuity: p(V • VC a ) = (2.4) 

State: p * pRT £ (C /W ) (2.5) 

a 

where the tc<tal enthalpy per unit mass is given by 

“ T - H+v2 /2 ( 2 . 6 ) 


In the above equations, * V*q R is the net rate of radiant energy absorbed 
per unit volume per unit time, K a represents the net rate of production of 

species a per unit volume per unit time, and W is the molecular weight of 

a 

species a . 

As a result of increased fluid enthalpy, the entire flow field in the 
precursor region is perturbed. By following the small perturbation technique 
of classical aerodynamics, the flow properties can be expressed in perturbation 
serie* as [11,12,22,23] 
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SHOCK LAYER 



Figure 2.1. Physical model and coordinate system 
for a Jovian entry body. 



P 00 (l + pj^ + P 2 + . . .) 

(2.7a) 

PooU + Pi + ?2 + • * 

(2.7b) 

v.(k + v x + v 2 + ...) 

(2.7c) 

H . + v »< H l + H 2 + •••> 

(2 . 7d) 

T . + T l + T 2 + •" 

(2.7e) 

■ c o +c a, + C cc, + ••• 

(2.7f) 


In these equations, all the perturbation variables (except temperature) 
are expressed in nondimensional form. The unit vector ic represents the 
direction of unperturbed free-stream velocity. 

If (L and K can be considered as first-order perturbation terms, 
then substitution of Eqs. (2.7) into Eqs. (2.1) -(2. 6) results in the first- 
order perturbation equations as 


Continuity: 

V • V x + 3p 1 /3z * 0 

(2.8) 

Momentum: 

3V x /3z « -(l/yM*)Vp 

(2.9) 

Energy: 

3 H T i /9z * V (p » V »> 

(2.10) 

Species: 

" V< P co V 

(2.11) 


H T. “ H 1 + V lz 

(2.12) 


The boundary conditions are that perturbation quantities vanish at z 
and that no singularities exist except at the origin. 


10 


It can be shown that the flow under consideration is irrotational [11,22]. 
Thus, there exists a potential 4> such that 

V 1 = V4> (2.13) 

For z-direction, integration of Eq. (2.9) results in 

P x * -(^) 34>/3z = -(YM4 )V 1z (2.14) 

Equation (2.8) now can be expressed as 

V 2 <J> + Sp^/Sz * 0 (2.15) 

In order to evaluate 3p^/3z and to relate to other variables, it is 

necessary to consider the gas model and radiation. 

For the Jupiter's atmosphere, the gas model is taken to be: C u * 89% 

H 2 

and Cjj * 11% by mole fraction (or ■ 80.82% and “ 19.18% by mass 

fraction) . The radiation effect on the gas ahead of the shock produces H+ , 

H, and electrons e“ by photodissociation and photoionization, and also 
increases the enthalpy. Any other species which may be produced are neglected. 
The contribution of radiation to the gas pressure is neglected. It is further 
assumed that the internal degrees of freedom of various species (i.e., vibra- 
tional and electronic modes) are not excited. For this gas model, the equation 
of state (for the first order perturbation) can be expressed as 

Pl - (400/180.17) [(C R + C H ^ + )/2j + (Tj/Tj + P x (2.16) 

By following the procedure described by Smith [11,12], the first-order 
perturbation relation for enthalpy is found to be 

H x * U/V*) 1 1 * 52 7 RT X + [(5/4)1^ + I/2]c h ^ + 

+ [(3/4)RT x + d]c h | (2.17) 
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where I and D represent the Ionization and dissociation energy respectively. 
It should be pointed out here that D in the above equation actually represents 
half the energy required for dissociation. 

As pointed out earlier, the upstream gas absorbs the energy radiated from 
the shock layer in the ultraviolet continuum range. The radiation from the 
perturbed gas due to recombination (i.e., emission) is neglected. The amount 
of radiative energy absorbed by the perturbed gas per unit volume and time, 

Q r , is given by 

X OO 

.. 2 w H v a(v)dv (2.18) 

where is the number density of H 2 , H v is specific irradiance and a(v) 

is the photon absorption cross section of H 2 at frequency v . 

In determining the rate of production of species in the precursor region, 
only photodissociation and photoionization are considered. Recombination is 
assumed to be a second-order effect and, therefore, is neglected in the present 
linearized treatment. The net rate of production of species, therefore, is 
given by [11,24] 


*H 


Q 1 N H X” (H v /hv) a D< v > dv 
2 


h 


a i 



(H v /hv) Oj(v)dv 


(2.19a) 


(2.19b) 


where m^ represents the weight of a H 2 molecule (in grams per molecule), 
and CJp(v) and <jj(v) are the absorption cross section for photodissociation 
and photoionization, respectively. 

Since the problem treated here is linear, it is permissible to obtain a 
solution for arbitrary frequency, and then integrate this solution over tne 
spectrum to obtain the solution for the problem. Thus, in the development that 
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follows, flow-field perturbations will be considered for a unit frequency 
interval. Equations (2.10) and (2.11) now can be written as 


3H„ / 3z - N„ 0(v)/(p V^) H 

*1 2 

(2.20) 

K h /3z ■ [.j » Y d 0(v)/(p. V. 

(2.21a) 

sc h 2+ /31 ■[ m i s h 2 y i "('''/(I 3 . H h v 

(2.21b) 


where Y^ and Yj represent photodissociation and photoionization yields, 
respectively. 

In order to express the governing equations in terms of perturbation 
potential, p^ is first eliminated by combining Eqs. (2.14) and (2.16). The 
resulting equation is then differentiated with respect to z and use is made 
of Eqs. (2.21). Next Eqs. (2. 12, 2. 13, 2. 17, and 2.20) are combined to give 

3p 1 /3z - -r 3 2 <J>/3z 2 - P v H v (2.22) 

where 

r - 0.727 y M 2 (2.23a) 

00 

p v • a v + b v /hv (2.23b) 

a v “ N H 2 a(v)/(p « V « H «> < 2 ‘ 24a > 

b v - -(a y m 1 /2)[(I - 0.89 RTJYj + (2D - 1.89 RT OT )Y D ] (2.24b) 


Upon combining Eqs. (2.15) and (2.22), the governing equation for the flow is 
obtained as 


v* - r 3 2 $/3z 2 - P H 

» y j v 


(2.25) 
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For axisymmetric case, this is expressed as 


r” 1 9/3r (r3*/3r) - T 3 2 4>/3z 2 - P v (2.26) 

Equations (2.25) and (2.26) are seen to be the classical potential 
equations for compressible flow with a forcing term proportional to radiation 
added. It should be pointed out that the form of Eq. (2.25) and (2.26) 
will be retained for any linearized gas model, although the expression for 
P v will depend on the gas model used. The potential for the flow induced 
by a radiant source with a spectral distribution is obtained by integrating 
the contributions of each frequency as 

$ «f " 0 dv (2.27) 

o 
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3. SOLUTIONS FOR SPECIAL CASES 


As discussed by Smith [11,12], solutions of the governing equations, 
presented in the previous section, can be obtained in special cases depending 
on the model used for the distribution of spatial radiation. If the radius of 
the radiating gas cap, R c , is large compared to the photon mean free path, 
then the problem can be treated like radiation from a plane source. On the 
other hand, when the radius of the radiating gas cap is small, then the 
problem can be treated like a spherical point source for radiation from the 
gas cap and a cylindrical point source for radiation from the wake. Note that, 
in general, R c may not be the same as the radius of the bow shock, R g . 

3.1 Radiation from a Plane Source 

For radiation from a plane source, it is essential to integrate the 
contribution over the plane, as attenuated by passage through the absorbing 
medium. The relation for H v , in this case, is given by [24] 

H V * 2 V 0) E 2 ( " k v z) (3.1) 

where q v (0) is the spectral radiative flux density at the shock wave, < v 

is the spectral absorption coefficient, and E n (t) is the exponential integral 
of order n . The expression for < v (which may also be interpreted as 
inverse of the photon mean free path) is given by 

< v - Njj O(v) (3.2) 

In this form < v represents the absorption coefficient of H~ molecules. If 
the number density Nj^ (and hence <^) can be taken to be Independent of z 
(which is a good approximation for small ionization and dissociation), then the 
optical depth is defined by 

S " V (3.3) 
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For the plane radiating source (where ^x,y $ * 0), therefore, a combination 
of Eqs. (2.25,3.1-3.3) results in 


-3 2 <J)/3C 2 - [2 P v q v (0)/(r k 2 )]E 2 (-C) 


(3.4) 


Integration of this yields the result 


<P » -12 P v q v (0)/(r< 2 )]E /( (-O 


(3.5) 


where the boundary condition of (9<J>/3?) -*■ 0 as C - 06 has been used. 

From Eq. (2.13), the velocity perturbations, ahead of the shock front, now 
can be written as 


v. » v, « 0 
lx ly 


(3.6) 


V 1Z - -< 2 P V iv«»/(rK v )iE 3 <-c> 


(3.7) 


From Eq. (2.14), the expression for pressure perturbation is found to be 


Pl " 12 Y P v q v (0)/K v ] (:^/DE 3 (-?) 


(3.8a) 


For high speed entry, M 2 » 1 and (M 2 /D 1 . Thus, Eq. (3.8a) can be 

00 oo — 

approximated by 


Pl - t2 Y P v q v (0)/< v ]E 3 (-O 


(3.8b) 


The expression for density perturbation is obtained by combining Eqs. (2.15) and 
(3.4) and integrating the resulting expression such that 


P 1 ■ 12 ? v V 0)/(I V !2 3 ; -« 


(3.9) 


By combining Eqs. (2.20,3.1-3.3), one obtains 


- [2 q v (0)/(P ao V£)]E 3 (-C) 
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Integration of this over the plane of the shock wave gives the result for the 
total enthalpy as 

H ? 1 " l 2 V 0)/(p » V »)3E 3 (-C) (3.10) 

By employing Eqs. (2.12,3.7, and 3.10), the expression for the static enthalpy 
is found to be 

H 1 " 2q v ( °){fl /(p 00 Vi)] + P v /( r < v )} E 3 (-0 (3.11) 

The concentrations of H and are given by integration of Eqs. (2.21) as 

C H “ l 2 w h 2 V (p « V » hv) l Y D (v) < J v ( 0)E 3 (-C) (3.12a) 

C H 2 + * [2 \ °1 /(P « V » hv)]Y x (v) q v (0)E 3 (-O (3.12b) 

By employing Eqs. (3. 8, 3. 9, and 3.12), Eq. (2.16) is solved for the temperature 
variation. For this case now all the flow properties at any point upstream of 
the shock can be determined. 

3.2 Radiation from Spherical and Cylindrical Point Sources 

The physical model for radiation from spherical and cylindrical point 
sources is shown in Fig. 3.1. A spherical point source is a source which 
radiates equally in all directions. A cylindrical point source is a source 
which radiates as a cylinder of infinitesimal radius and length. For both 
cases, the incident radiation at any field point s is given by [11,12] 

H v " (A \/ s2 ) ex P( _< v s ) ( sin (3.13) 

In this equation, A^ represents the radiative strength of the source, s is 
the distance from the source, and 9 is the angle between the free stream 
velocity vector and a line from the field point to the center of the source. 
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Figure 3*1. Approximation of radiating flow by point sources. 



The superscript j ■ 0 for a spherical point source and 1 for a cylin- 
drical point source. 

Equation (3.13) can be substituted In Eqs. (2.25) and (2.26) to obtain 
the corresponding equations for the perturbation potential. Within the con- 
fines of the assumptions made In obtaining Eq. (3.13), however, both problems 
(spherical as well as cylindrical point source) can be considered to be 
axisymmetrlc . The governing equation for the perturbation potential, there- 
fore, can be written as 

rf 1 3/3n (n3<J>73n) - r 3V/3? 2 - A U' 2 e‘ u (sin 6) J (3.14) 

where 

U - < v s, n * < v r, - k v 4>, A - A(v) - < v P v A v 

A procedure for general solution of this equation is suggested by 

Smith [11]. For entry flows, however, M 2 >> 1 and Eq. (3.14) can be solved 

00 

by expanding <p' in a series in (1/0 in the vicinity of the body. Thus, 
one can express <p' as 

<t>' - -(< v /r)[Fj(;,n) + d/r)F j (1) (c,n) + (I/O 2 Fj (2) (c,n) + ...J (3.15) 

where Fj’s are functions for perturbation potential. Substitution of this 
relation into Eq. (3.14) gives 

3 2 Fj/ 3C 2 - U~ 2 exp(-u) (sin '?) j (3.16) 

and 

3 2 F j (n) /3C 2 - -rf 1 3/3n(n 3Fj (n " 1) /3r l ) (3.17) 

The problem, therefore, is reduced to quadratures in the vicinity of the 
body. In the present analysis only the term in (1/D will be retained. By 
integrating Eq. (3.16) twice the expression for F^ is obtained as 

F.(C,n) - /°U~ 2 exp(-U o )(n/U o ) J (’ - ; o )dC o (3.18) 
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where U« • n 2 + ? 2 . For convenience, let us denote 

O 

Gjtt.n) • SF.j/3; •f' u -2 exp(-y o )(n/U o ) J d? £ 


(3.19) 


^(c.n) - 3F./3n - u" 3 exp(-u o ) |(n/u)^ tn + 

** * —00 

(n/y 0 )(2 + j)l - j)(? - ; 0 )d; o 

With these definitions of F j > > and H j > the perturbation quantities 

can be expressed as 


(3.20) 


4>' - -(A/OFjU.n) 


(3.21) 


v lr * (A/r)Hj(C,n) 


(3.22) 


v lz - -(A/DGjU.n) 


(3.23) 


P x - yAGj(;,n) 


(3.24) 


0 1 - (A/DG j (;,n) 


(3.25) 


^ “ (K v V p » v*)Gj(;,n) (3.26) 

C H " (m l A v K J /p » V oo hv)Y D (v) Gj(5,n) (3.27a) 

C H 2 + " ( ®1 \> V « hv)Y i (v) G^r.n) (3.27b) 

Note that for the case of spherically radiating point source j » 1 in the 
above equations. Also, these equations are obtained for arbitrary frequency. 

The expression for total potential, for this case, can be obtained by combining 
Eqs. (2.27) and (3.15) as 


$ - ~(l/n/ 0 “ [AOO/KjFjC.rOdv 


(3.28) 
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Furthermore, it should be noted that the above solutions are valid in the 
region where [y“ 2 exp(-y) (sin 9)^] does not vanish. This is the case 
of spherically symmetric flow ahead of the entry body and is of primary 
concern in the present study. Other cases involving cylindrical point source 
are discussed in [11,12]. 

3.3 Solutions for the Transition Range (^Rc " 0(D) 

As mentioned earlier, for <v R c << 1 * the axisymmetric solutions for 
the spherically and cylindrically radiating point source are valid and for 
< V R C » 1 , the one-dimensional equations apply. The solutions for the spheri- 
cally radiating point source approach one-dimensional solutions as < V R C 00 • 
Thus, spherically radiating point source solutions are valid for the precursor 
flow ahead of a blunt body with < V R C >:> * an£ * K v R c <<: * * Since the 
< V R C » 0(1) range lies between these two limits, the spherically radiating 
point source analysis could be applied (to a good approximation) also in 
the transitional range. 


4. PERTURBATION EQUATIONS FOR PHOTOABSORPTION MODEL 


In order to obtain specific results for the perturbed quantities, it 
is essential to have a realistic model for the spectral absorption coefficient 
of the absorbing gas. The photoabsorption model employed in this study is 
discussed in this section and the governing perturbation equations are expressed 
in special form for this model. 

4.1 Photoabsorption Model 

Photoionization absorption coefficient is a continuous non-zero function 
of photon energy (because of bound-free transition) for all values of photon 
energy that exceed the ionization potential of the atom. Similar remarks 
apply to the photodissociation and radiative recombination. A critical 
review of ultraviolet photoabsorption cross secticis for molecules of astro- 
physical and aeronomic interest, available in the literature up to I.'-?l, is 
given by Hudson [25] . Specific information on photoionization and absorption 
coefficient of molecular hydrogen is available in [25-31] . 

Photoionization and absorption cross sections of H£ , as obtained from 
references [25-31], are plotted in Fig. 4.1. From this figure it is evident 
that the ionization continuum starts at about 804 A and continues towards 

e o 

lower wave lengths. Between the wave lengths of 600 A and 804 A, the absorp- 
tion cross section for ionization continuum are included in the total absorp- 
tion (i.e., absorption due to ionization as well as dissociation). For wave- 

O 

lengths below 600 A, however, the ionization continuum absorption is equal to 
the total absorption. The total absorption cross section for continuum range 

P 

below 304 A can be closely approximated by the two rectangles (I and II) 
shown in the figure with broken lines. The ratio of ionization cross section 
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Figure 4.1. Absorption cross section of in ultraviolet region. 







to the total absorption cross-section (i.e., the value of Yj) is taken to 
be unity for rectangle I and 0.875 for rectangle II. For wavelength greater 

O 

than 804 A (where hv is below ionization energy) , the value of Yj is taken 
to be zero. Little information is available in the literature on the absorp- 
tion cross-section for dissociation of H 2 molecules. There is strong evi- 

Q 

dence, however, that photodissociation starts at about 2600 A and continues 

O 

towards lower wavelengths to about 750 A (26-28]. There are also a few 
diffuse bands in this spectral range [26,28]. Thus, it becomes difficult 
to evaluate the absorption cross section in this spectral range. For this 

O 

study, the absorption cross section in the spectral range between 804 A and 

O 

2600 A was approximated by the rectangle III. The specific values of a(v) 
for the three rectangles are found to be ^(v) - 4.1 E-18, ©^(v) - 8.2 E-18, 
and ^m(v) * 2.1 E-18. The value of Y^ is taken to be zero for rectangle I 
and 0.125 for rectangle II. 

4.2 Precursor Equations for the Photoabsorption Model 

In accordance with the photoabsorption model assumed in the previous 
subsection, quantities <J(v), A, Y j and Y^ are taken to be constant within a 
given frequency range. Thus, these quantities are assumed to have constant 
(but different) values for rectangles I, II, and III in Fig. 4.1. It is 
further assumed that the gas cap radiates as a gray body and the effective 
emissivity of the gas cap is given by 


e - q(0) / (O T**) 


(4.1) 


where q(0) is the radiative flux density from the shock layer and T g is the 
shock temperature. 

The spectral radiative flux at the shock front is given by the relation 


q v (°) - e TT B v (T s ) 


(4.2) 


2 ' 




In this equation, B..(T ) represents the Planck function which is given by 

v 3 

B v (T) * (2h/c 2 )(kT/h) 3 { v 3 /[exp(v) - 1]} (4.3) 

where v ■ hv/kT . 

It now remains to obtain the relation for in Eq. (3.13). By noting 

that at the shock front » q v (0) and s » R^ » a rear r an 8 eoent of 
Eq. (3.13) gives 

\ - q v (0) R 2 exp(< v R c ) (4.4) 

With the above information, final relations for the total perturbation 
quantities now can be obtained. For spherically radiating case, the perturba- 
tion quantities can be expressed by an equation of the form [11,12] 

= Qf* A(v) Gj(4,n)dv (4.5) 

where ip represents V lz * Pi , or , and 8 represents the factor ahead 
of the integral. The radial component of the velocity perturbation takes 
this same form but Gj is replaced by Hj . By employing the definition of 
A(v) , and combining Eqs. (2.24) and (4.5) results in 

'P x SR \f Q < V U V + (b v /hv) ] exp(K v * c ) q y (o) cv (4.6) 

Since a v , b v , and are assumed to be constant over sections of the 

wavelength range, Eq. (4.6) can be expressed as 
n 

’P * BR 2 £ k. exp(K. R ) G (k.z, < r) * 

<- i j. c j i 1 

( a i/ V2i q v (0)dv + (b /h) f^ 2i [q v (0)/v]dv| (4,?) 

V li V li 

The two Integrals in the above equation are evaluated by using Eqs. (4.2) and 
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- > 


(A. 3). The first one can be written as 


/ V V 

2i q (0)dv = e TT(2h/c 2 )(kT /h) 4 f 2i [v 3 /(e V - l)]dv 

, , V S 

V li v li 

By employing the definition of a and denoting 

N^(v) -X* [v^/(e V - l)]dv 

Eq. (4.8) can be expressed as 

r v 2i 

J q v (0)dv - (15/ir 4 ) e a T 4 [N^v.^) - N 3 (v 21 )] 

V li 

Similarly, the second integral in Eq. (4.7) can be expressed as 
r v 2i 

J [q v (0)/v]dv = (15/7T 4 ) EOT 4 (h/kT ) [N- (v. . ) - N,(v 0 .) ] 

\) S S Z 11 Z ZX 

v i j 


(4.8) 


(4.9) 


(4.10) 


(4.11) 


Substitution of Eqs. (4.10 and (4.11) into Eq. (4.7) gives the final relation 
for ip from which the perturbation quantities can be determined. The 
quantities lij , , 0^^+ have the form of Eq. (4.7) with the which 

appears ahead of the exp(<jR c ) squared. 

In order to write specific relations for the perturbation quantities, 
it would be convenient to define 


8 X - ( 15 /tt 4 ) [q(0)/r] 

(4.12a) 

B 2 - (15/tt 4 ) y q(0) 

(4.12b) 

S 3 « (15/tt 4 ) [ q ( 0) / ( Poo V^)] 

(4.12c) 

8 4 - (15/tt 4 ) [q (0) m 1 /(p oo VJ J 

(4.12d) 

I(vJ) - f 21 | v p / [exp(v) - l]|dv 
li ’ 

(4.13) 

B(a i ,b i ) - a ± I(v 3 ) + (b^kTj I(vJ) 

(4.14) 
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For the spherically radiating case, the perturbation quantities now can 
be expressed as 




n 

S * S 4 R c £ Y D K i exp(K i R c ) G j ( ^> I(vJ) (4. 


\+ - e 4 R* Z ^ MPCK, R c ) GjCc.n) I(vj) (4. 

T 1 " T » I Pi * p l “ (200/180.17) (C H + C H )] - (4. 

For the plane radiating source, the perturbation quantities can be 
written as 


v a ‘ - 2S i £ ''I 1 E 3 ( -i> B <V b i> < 4 - 

H 1 ’ 26 3 £ E 3<-5.) (4. 

J 1-1 J 1 

C H ■ 28 4 £ \ E 3<- C i ) < kT s>" I(v l> «• 

C H 2 + ■ 28 4 t \ V'V (k V‘ 1 < v! i> (*• 


The expressions for and in this case are the same as for the 

spherically radiating source but care should be taken In using the right 
relations for V^ 2 , 0^+ and . 

Depending on the order of optical thickness, either spherical or plane 
radiating source relations are employed in actual calculations of the per- 


turbation quantities. 
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5. DATA SOURCE AND SOLUTION PROCEDURE 


Jupiter's atmospheric conditions, as obtained from reference [32], 
are shown in Fig. 5.1 for different altitudes. The temperature of the atmo- 
sphere (i.e., T m ) is taken to be constant at 145°K and the enthalpy is given 

by H_ ■ 1.527 R T . The entry velocity range considered are between 

00 

28-45 km/sec. The value of R c * 25 cm was assumed for this study. The 
number density of H 2 can be computed by using the relation 


N 


H 


2 


(7.2431172 x IQ 22 ) (PjTjX, 

2 


(5.1) 


where Xjj 2 is the mole fraction of H 2 and has units of N/m 2 . 

Free stream and shock conditions used in this study are listed in 
Table 5.1. Shock temperature and q(0) -values were calculated by employing 
the computational procedure developed by Sutton [4] and Moss [9]. Since 
viscous effects are pronounced primarily in the vicinity of the body, only 
inviscid shock layer formulations were considered in calculating T g and 
q(0). Further details on shock layer solutions and effects of shock precursor 
heating on radiative flux to the body are given in a separate report [33]. 

Before evaluating actual values of the perturbation quantities, it was 
considered necessary to investigate the range of different intervening para- 
meters of the governing equations. The values of the inverse mean free path 
(i.e., < v » a(v) N|j 2 ) and the product < v R c were calculated for the Jovian 
entry conditions and these are shown in Fig. 5.2 for the three different 
values of the photoabsorption cross section of Fig. 4.1. From this figure it 
is evident that the product (< R c ) >> 1 in most cases of interest for the 
Jovian entry. Thus, one could employ only the plane radiating source formula- 
tions for determining the perturbation quantities in the precursor region. 
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Table 5.1 Free- stream and shock conditions 
for Jovian entry. 


Free stream 

V , km/sec 
00 

, 

T , °K 
8 

q(0) , erg/cm 2 

Z ■ 95, km 

38 

16,610 

■ 

1.35 E12 

p« - 1.29 E-3, kg/m 3 

35 

15,400 

7.75 Ell 

P m - 673, N/o 2 

32 

14,080 

3.52 Ell 


30 

13,550 

2.01 Ell 

Z - 103 

40 

16,890 

1.16 Ell 

Po, - 8.56 E-4 

35 

15,040 

4.70 Ell 

P» - 448 

33 

14,250 

3.28 Ell 


30 

12,810 

1.142 Ell 

: 

2 - 116 

45 

18,227 

1.09 E12 

p m - 4.65 E-4 

39.09 

15,886 

4.76 Ell 

Poo “ 244 

35 

14,480 

2.18 Ell 


30 

12,480 

4.87 E10 

Z “ 131 

43.21 

16,390 

3.86 Ell 

- 2.32 E-4 

38 

15,210 

1.61 Ell 

Pc, - 122 

35 

13,880 

8.72 E10 


30 

12,030 

1.90 E10 

Z - 150 

42 

15,050 

9.60 E10 

) 

Pgo - 9.29 E-5 

40 

14,520 

6.96 E10 

Poo - 49 

35 

13,140 

2.57 E10 


30 

1 

11,600 

I 

6.20 E9 























INVERSE MEAN FREE PATH (A), cm' 1 



Figure 5.2. Inverse mean free path and (21^) as a function of altitude for the 
three regions of photoabsorption model. 
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Values of the intermediate functions a v and (defined by Eq. (2. 24)) 

were calculated for the photoabsorption model at different altitudes. Since 
To,, is taken to be constant, the values of functions (a^ V^) and (b v V,,,) 
were found to be constant for all altitudes. Another intermediate function, 
A(v), is defined in Eqs. (3.14) and (4.5). For < v R » 1 , the expression 
for A(v) can be written as 

A(v t ) - (15 /tt 4 ) q(0)/< vi E 3 (0) B(a i ,b i ) (5.2) 

The values of A(v^) were calculated for the three spectral range of the 
photoabsorption model and the results are illustrated in Fig. 5.3. For any 
altitude, the value of this function increases with increasing entry velocity. 

The set of Eqs. (4.15)-(4.22) for spherically radiating source and 
Eqs. (4.23)-(4.26) for plane radiating source can be solved numerically to 
obtain the perturbation quantities. As mentioned earlier, for Jovian entry 
conditions, it is necessary to solve only the set of equations belonging to 
the plane radiating source case. 

For multicomponent systems, it is physically realistic (and a general 
practice) to define the total enthalpy of the gas entering the shock layer 
by the relation 

Hj - H q + V 2 /2 (5.3) 

In this equation, V represents the local fluid velocity, and H Q is 
referred to as the absolute enthalpy and is equal to the sum of sensible 
enthalpy and chemical energy at 0 8 K [34,35]. In terms of the perturbation 
velocity, the local fluid velocity is given by 

V - V (1 T V. ) (5.4) 

For a multicomponent system, the expression for H 0 is given by 

H - T C II (5.5) 

o *-* a o r 
a a 
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Figure 5.3b. A(v) as a function of free-stream velocity for 
spectral range II. 
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where c^ is the mass fraction of species a . For any species a , the 

relation for H is given by [34,35] 

°a 

H « RT[a 1 + (a 2 /2)T + (a 3 /3)T 2 + (a 4 /4)T 3 
°a 

+ (a 5 /5)T 4 + a 6 /T] (5.6) 

where R is the universal gas constant (» 1.98726 cal/mole - °K) and 
T is the local fluid temperature in °K. For different species, values of 
constants a^, a 2 >...ag are given in [35], and for species under present 
investigation they are listed in Table 5.2. 

It should be noted that values of Hj can be calculated by considering 
or neglecting the precursor effects. When precursor effects are considered, 
then Hj, is defined by 

»T ' < H t> PE " “VPE + V2/Z (5.7) 

For the case with no precursor effects, H^. is given by 

lI T “ (1 V«PE = * H o*NPE + V » /2 < 5 ' 8 ) 

It should be emphasized 'ere that, for the case with no precursor effects, 
the temperature in Eq. (5.6) is the free-stream temperature T^ . 

The per cent difference in the total enthalpy with and without the 
precursor effects can be expressed by 

% - PD - |[(H T ) pE - <H T ) NpE l/(H T ) NpE } x 100 (5.9) 

It should be noted that an appropriate value of H^. is needed to determine 
the conditions inside the bow shock by employing the Rankine-Hugoniot relations. 
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Table 5.2. Temperature coefficients for thermodynamic functions for hydrogen species. 


Substance 


H 


2 


H 


+ 

2 


H 


Temperature 
Interval, °K 

Coefficients 

a i 

a 2 

a 3 

a 4 

a 5 

a 6 

1000 < T < 5000 

3.0436897 

6. 1187110E-4 

-7 . 3993551E-9 

-2.0331907E-11 

2.459379E-15 

-8.5491002E2 

300 < T < 1000 

2.8460849 

4.1932116E-3 

-9.6119332E-6 

9.5122662E-9 

-3.309342E-12 

-9.6725372E2 

1000 < T < 5000 

3.3287156 

2.5050678E-4 

1.4224521E-7 

-4.4590247E-11 

3.733756E-15 

1.7997470E5 

300 < T < 1000 

2.817375 

3. 657610E-3 

-7.9655480E-6 

8.26J4000E-9 

-3.090228E-12 

1.8002 739E5 

1000 < T <5000 

2.500 

0 

0 

0 

0 

1 2.5470497E4 

300 < T < 1000 

2.500 

0 

0 

! 

0 

0 

2.5470497E4 



6. PERTURBATION RESULTS 


The flow perturbation quantities V lz * Pi j j Cjj t i T^ ) and 

were calculated numerically and the results are illustrated in Figs. 6.1 - 

6.12. In Figs. 6.1 - 6.6, perturbation quantities are shown as a function of 

distance from the shock for different altitudes and a constant entry velocity 

of 35 km/ sec. The first set of curves (Figs. 6.1a - 6.6a) are plotted against 

the nondimenslonal distance z/R , while the second set (Figs. 6. lb-6. 6b) are 

c 

plotted as a function of the physical distance z from the shock. In Figs. 

6.7 - 6.12, the perturbation quantities (just ahead of the bow shock) are 
illustrated as a function of the free-stream velocities. Since » 

separate results were not illustrated for the density perturbation. From these 
figures it is evident that the magnitude of perturbation quantities, in general, 
depend on the distance from the shock, altitude of entry, and entry speeds. 

Figures 6.1 - 6.6 show that at a fixed entry velocity, the perturbation 
effects are greater for lower altitudes and at locations just ahead of the 
shock. This, however, would be expected because the number densities of parti- 
cipating species are greater at lower altitudes and at these altitudes most 
radiative energy from the shock gets absorbed in the immediate vicinity of the 
shock front. At higher altitudes, perturbation effects are significant to a 
larger distance from the shock front. This is because, at these altitudes, the 
number densities of participating species are small and radiation effects are 
felt farther into the free-stream. 

Specific results presented in Figs. 6.1b - 6.6b indicate that use of the 
small perturbation theory is justified in determining the velocity, density, 
mass fraction and total enthalpy variations. These variations are small because 
at high entry speeds, the gas has not had enough time for expanding. 
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DISTANCE FROM SHOCK (r/R c ) 


Figure 6.1a. Velocity perturbation as a function of distance from 
the shock at different altitudes and a constant free 
stream velocity. 
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Figure 6.2a. Pressure perturbation as a function of distance from 
the shock at different altitudes and a constant free- 
stream velocity. 
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Pressure perturbation as a funct 
the shock at different altitudes 
stream velocity. 








DISTANCE FROM SHOCK (z/Rq) 


Figure 6.3a. Mass fraction of H as a function of distance from 

the shock at different altitudes and a constant free- 
stream velocity. 
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OISTANCE FROM SHOCK, cm 

Figure 6.3b. Mass fraction of H as a function of distance from 

the shock at different altitudes and a constant free- 
stream velocity. 
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DISTANCE FROM SHOCK (z/R c ) 


Figure 6.4a. Mass fraction of H^+ as a function of distance from 
the shock at different altitudes and a constant free- 
stream velocity. 
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For example, just ahead of the shock, the value of (V/V^) is 0.9992 for 

Z ■ 95 km and is equal to 0.99975 for Z ■ 150 km. Similarly, * 

6.8 x 10“ 3 for Z - 95 km and - 2.4 x 10~ 3 for Z - 150 km (i.e., 

0.68% increase in total enthalpy at 95 km and 0.24% increase at 150 km). The 
static pressure and temperature variations, however, cannot be considered small 
This is because for 2 * 95 km, p^ ■ 2 and T^ * 300°K, and for Z * 150 km, 

Pj * 0.64 and T^ * 94 °K . For these variations, therefore, one could ques- 
tion the validity of the small perturbation theory. 

For different altitudes of entry, perturbation results (just ahead of the 
shock) are illustrated in Figs. 6.7 - 6.12 as a function of entry velocities. 
The results are shown only for the range of entry velocities for which free- 
stream and shock conditions are available (see Table 5.1). These results again 
indicate that the perturbation effects are greater for lower altitudes. As 
would be expected, for any specific altitude, the effects are larger for higher 

entry velocities. This is a direct consequence of greater radiative energy 

transfer from the shock to the free-stream at high entry speeds. For the most 
part, variations in the velocity, mass fractions, and total enthalpy again are 
seen to be small. For example, for an entry body at an altitude of 95 km, 
the total enthalpy of the gas (H^) entering the shock wave is increased from 
about 0.68% at V * 35 km/sec to 1% at V ■* 38 km/sec. For Z * 150 km, 
however, increases from 0.24% at 35 km/sec to 0.66% at 42 km/sec. The 

variations in the static pressure and temperature, in some cases, are seen to 
be several times greater than the ambient values. These large variations, 
however, occur for conditions where dissociation is high and the validity of 
the entire theory is questionable [11,12]. This point is discussed further in 
the next section. 
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Figure 6.7, Velocity perturt 
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(just ahead of the shock) as a function of free-stream 
altitudes . 
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:urbation (just ahead of the shock) as a function of free “Stream 
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The results of absolute enthalpy (as defined by Eqs. (5.5) and (5.7)) are 
plotted in Fig. 6.13 as a function of free-stream velocities. The trend of 
results of this figure are essentially the same as for specific total enthalpy 
perturbation results shown In Fig. 6.12. This again illustrates that the 
Increase in enthalpy due to precursor absorption are greater for lower alti- 
tudes and higher entry velocities. 

The per cent difference in total enthalpy of the gas (with and without 
precursor effects) entering the shock wave is illustrated in Fig. 6.14 for 
different altitudes. The results indicate that the maximum increase in total 
enthalpy is about 1.6% for Z « 95 km and entry speeds of 38 km/ sec. For other 
entry conditions, the changes are seen to be smaller. 

A few conclusions can be drawn from the results presented in this section. 
Within the limitations of the small perturbation technique used in this study, 
the results indicate that variations in velocity, density, species concentra- 
tion, and enthalpy are small as compared to perturbations in pressure and temp- 
erature. Precursor effects, in general, are greater for lower altitudes and 
higher entry velocities. At higher altitudes, however, precursor effects are 
felt farther in the free-stream. At any particular altitude, the effects 
Increase with increasing entry velocities. Specific results Indicate ‘hat for 
Jovian entry velocities lower than 28 km/sec and altitudes of entry higher than 
95 km, the precursor effects definitely can be neglected. For other entry con- 
ditions, the extent of flow perturbations and its influence on the entire flow 
field ahead of the entry body should be investigated. 
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Figure 6.13. Absolute enthalpy (just ahead of the shock) as a function 
of free-stream velocity for constant altitudes. 
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FREE-STREAM VELOCITY (V„), km/sec 

Figure 6.14. Increase in total enthalpy (just ahead of the shock) as a function of 
free-stream velocity for constant altitudes. 








7. ALTERNATE APPROACH: THIN-LAYER APPROXIMATION 





Another approach to investigate the precursor effects will be to apply 
the concept of thin shock layer theory (usually applied to hypersonic shock 
layer flows (36,37]) to the precursor region flow field. For this purpose, 
a curvilinear orthogonal coordinate system, shown in Fig. 7.1, is selected. 
In thi 3 figure s is the distance (measured from the stagnation point) 
along a reference surface (body or shock) and n the distance along the 
normal to this surface. For convenience, the reference surface is taken to 
be the outer edge of the shock layer. 

The differential equations for a hypersonic plane or axisymmetric flow 


can be written in the present coordinate system as [36] 

(3/3s)(pur*) + (3/3n) (pvXr^ ) * 0 (7.1a) 

p(u(3u/3s) + Xv(3u/3n) - Kuv] + (3p/3s) * 0 (7.1b) 

p(u(3v/3s) + Xv(3v/3n) + Ku 2 ] + X(3p/3n) » 0 (7.1c) 

p((u/X)(3H/3s) + v(3H/3n) ] + (Xr J )’ 1 [ (3/3n) (Xr J q_) ] - 0 (7. Id) 

P((u/X) (3C a /3s) + v(3C a /3n) - ■ 0 (7.1e) 


where K - K(s) - 1/R S , X » 1 + Kn , and j • o for plane flows and 
1 for axisymmetric flows. 

If the precursor region is assumed thin, then one can make the approxi- 
mations that (n/R g ) << 1 , 3/3s << 3/3n , and is not a function of n . 

In this case X • 1 , and Eqs. (7.1) reduce to 
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y 





Figure 7.1. Curvilinear orthogonal coordinate systems for 
thin-layer approximation. 
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(3/3n)(pv) - 0 


(7.2a) 


pv(3u/3n) ■ 0 

(7.2b) 

pv(3v/3n) + (3p/3n) * 0 

(7.2c) 

pv(3H/3n) + (3q^/3n) * 0 

(7. 2d) 

pv(3C a /3n) - K 0 - 0 

(7.2e) 

Direct integration of these results in 


P V - P. V oo 

(7.3a) 

p„ v a, <“ - ■ 0 

(7.3b) 

P„ <* ■ »„> + (p - p„) ■ 0 

(7.3c) 

p. ». CH - HJ ♦ <!„ • o 

(7.3d) 

»„ »„ - K 0 - 0 

(7.3e) 

where it has been assumed that q R - o . 

OO 

In present application to the hydrogen-helium atmosphere. 

Eq. (7 . 3e) 

will be written for atomic hydrogen and hydrogen ions. In Eq. 

(7.3d), H 

represents the total enthalpy and is given by the relation 


H ■ H t - h + (u 2 + v 2 )/2 

(7.4) 
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where 


h - 1.527 RT + [ (5/4)RT + I/2]C„ ^ 

h 2 + 

+ [(3/4)RT + D]C h (7.5) 

Note that Eq. (7.5) is slightly different than the relation for perturba- 
tion enthalpy given by Eq. (2.17). 

For a diffuse nonreflecting shock front, the expression for one-dimen- 
sional spectral radiative flux and its spatial derivative is given by [24,38] 

WV - 2 £ v *bv<V W + 2 X ' bv (t) E 2 (t v - c > dt < 7 - d > 

- d W dT v - “.S. (I s ) W - 2 V. (T) + 2 fo e bv (t) E 1 (T - t)d,; < 7 - 7 > 

where 

T v K v dn ' < 7 * 8 > 

o 

e bv (T) - tt B v (T) = 7T(2h/c 2 ){v 3 /[exp(h /kT) - 1]} (7.9) 

It should be noted that these equations do not account for any radiation 
from the free-stream. As before, if the number density of H2 is assumed 
constant, then in the above equations k v becomes independent of position. 

The expression for total radiative flux is given by 

q R (n) 'X” q Rv (T v )dV (7.10) 

For a gray shock front, a combination of Eqs. (7.6) and (7.10) results in 

q (o) " T 4 » which (as would be expected) is the same as q(o) defined 
R s 

in Eq. (4.1) . 

If emission from the cold gas in front of the shock is neglected, then 


for a gray shock front, one can write 




L 


Or<"> * 2 e / 0 " v <V 2 3 ( V d '' ‘ 2 /o” V 0) E 3 ( V' 


(7.11) 


-dq R /dn - 2 £ / o °° K v e bv (T s ) E^dV = 2 £ < y q v <0) E^dV (7.12) 


where q v (o) Is defined by Eq. (4.2). 

The general expression for the total radiative flux is obtained by 
combining Eqs. (7.6) and (7.10) as 

q R (n) “ 2 fo ( q v (0) E 3 (K v n) + ™\>fo B v (T) W n - n ')] d n'}dv (7.13) 

where is assumed to be independent of position. For the spectral model 

considered in section 4, Eq. (7.13) can be written as 

n v 

q R (n) - 27 T £ j( 15 / 7 T s ) q(0) E,(<.n ) f 21 [v 3 /(e V - 1) Jdv 

R i-1 1 31 v li 


+ K i-C E 2 [l< l (n - B v (I) dv dn 'l 


(7.14) 


where v * hv/k T c . The final form of the energy equation now can be 
s 

obtained through a combination of Eqs. (7.3d, 7.4, 7.5, and 7.14). 

Either by following the information given in Eqs. (2.18, 2.19, and 3.1) 
or from Ref. [19], the expression for concentration of the species a in 
Eq. (7.3e) can be written as 

P* v « (9c a /9n) * -®i f 0 °° [(3q Rv /9n)/hv]dv ( 7 .; 

The appropriate expression for (3qp^/9n) , in this case, is given by 
Eq. (7.12). Thus, Eq. (7.15) can be expressed as 

dC a /dn - (2m 1 /hp oo vjf* [< v q v <0) E « n)/v]dv (7.: 


By noting that < v * c(v) and following the procedure of section 4 
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the relations for individual species are found to be 

n 

. C H - -2 6 4 £ t„ t E 3 (T 1 )(kT s )-‘ KvJ) (7.17) 

\* - -2 e 4 t * 1± W^V 1 

where 8^ , I(v^ 2 ), and are defined by Eqs. (4.12d, 4.13, and 7.8) 
respectively. It should be noted that Eqs. (7.17) and (7.18) are exactly 
the same as Eqs. (4.25) and (4.26). 

By employing the governing equations presented in this section and the 
spectral information of section 4, numerical results were obtained for 
velocity, pressure and temperature variations for different values of n 
at s = o . Specific results for an altitude of Z » 116 km are compared 
in Figs. 7. 2-7. 4 with corresponding results of the small perturbation theory. 

As indicated earlier, the equations for species concentration in this case 
are found to be exactly the same as for the small perturbation case. For 
the range of parameters considered, the results for velocity, pressure and 
temperature obtained by the two procedures are seen to be in excellent agree- 
ment. It is obvious from these results that either approach could be utilized 
in the investigation of the precursor region flow field. It was noted in 
section 5 that for the Jupiter's entry conditions, the general governing 
equations of the small perturbation theory reduced to the case of simple 
plane source. As such, use of this method to Jupiter's entry case is restricted 
to one-dimensional analyses. The advantage of thin layer approximation pro- 
cedure is that it is physically more convincing and it can be extended easily 
to three-dimensional and axisymmetric cases. Furthermore, in more realistic 
situations, the thin layer approximation can be relaxed and the analysis can 
be extended easily to general cases. 
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8. CONCLUSIONS 


Governing equations have been presented for investigating the 
precursor region flow by employing the small perturbation theory of 
classical aerodynamics and thin-layer approximations of hypersonic flows. 

In small perturbation method, the perturbation velocity potential is 

found to be governed by the wave equation with a driving term due to radi- 

ation absorption, ionization and dissociation. The thin layer approxima- 
tion reduces the general hypersonic flow equations to simpler forms for 
which solutions are obtained by employing an iterative procedure. 

By employing appropriate thermodynamic and spectral data for the 
hydrogen-helium atmosphere, variations in precursor region flow quantities 
were calculated by the two entirely different methods. For Jovian entry 
conditions, one-dimensional results obtained by the two methods were 
found to be in good agreement for the range of parameters considered. The 

results, in general, indicate that for certain combinations of entry speeds 

and altitudes of entry, the precursor effects cannot be ignored while 
analyzing flows around Jovian entry bodies. Specifically, it is seen that 
at an altitude of 95 km, the precursor effects are important for entry 
velocities greater than 35 km/sec. 

The usefulness of the thin-layer approximation in analyzing the pre- 
cursor region flow is demonstrated. The main advantage of this method is 
that it is physically more convincing and its use can be extended easily 
to axisymmetric and three-dimensional cases. It is suggested that precursor 
region flow phenomena be investigated in general without making assumption 
of the thin layer approximation. It might even be advisable to modify the 
radiation model for the precursor region absorption. The extent of pre- 
cursor effects on the entire shock layer flow phenomena should be investi- 
gated thoroughly. It might even be essential to include two-dimensional 
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model for radiative flux and a detailed spectral model for radiation 
absorption and emission. 
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APPENDIX A1 


AO, A1 

A2, A3 

BO, B1 

B2, B3 

CH 

CHE 

CHH 

CH2 

DEN 

DEN 2 

ETH 

ETHP1 

ETHP2 

ETHP3 

ETHP4 

EX 

GARM 

GD 

GI 

GMARC 

GNH2 

P 

PRE 


EXPLANATION OF SYMBOLS USED IN THE COMPUTER PROGRAM "PERC" 


Exponential integral constants 


Exponential integral constants 

Mass fraction of H 
Mass fraction of H 

e 

Mass fraction of H 2 
Mass fraction of H 2 + 
Perturbation density 
Free-stream density, g/cm 3 


Sum 

of 

sensible 

and 

chemical 

enthalpy 

of 

mixture gas, cal/g 

Sum 

of 

sensible 

and 

chemical 

enthalpy 

of 

h 2 . 

cal/g 

Sum 

of 

sensible 

and 

chemical 

enthalpy 

of 

h 2 +. 

cal/ g 

Sum 

of 

sensible 

and 

chemical 

enthalpy 

of 

H, 

cal/g 

Sum 

of 

sensible 

and 

chemical 

enthalpy 

of 

H e’ 

cal/g 


Exponential integral 
Specific heat ratio 
Dissociation energy, erg/mole 
Ionization energy, erg/oole 
Mach number 

Number density of H, , molecule/cm 3 
Boltzmann constant, erg /°K 
Perturbation pressure 



PRI Free-stream pressure, dyne/ cm 3 

QO Radiative heat flurc at shock front, erg/cm 3 -sec 

RCAP Shock radius, cm 

ROE Absorption cross section, cm 2 

T Absorption coefficient, cm"* 

TF1 Perturbation temperature, °K 

TF111 Temperature, °K 

TH Planck's constant, erg-sec 

TO Free-stream temperature, *K 

TS Temperature at shock, *K 

U Universal gas constant, erg/mole-°K 

V Frequency, sec -1 

VINF Free-stream velocity, cm/sec 

VZ Perturbation velocity component normal to the shock surface 

W1 Molecular weight, g/molecule 

Z Distance from the shock 
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APPENDIX A2 


EXPLANATION OF SYMBOLS USED IN THE COMPUTER PROGRAM "THIN" 


AO, A1 
A2, A3 
ARF 
BO, B1 
B2, B3 
CH 
CHE 
CHH 
CH2 
CK 
DEN 
DENI 
DH 

El, E2, E3 

ETH, ETHi 

ETH2, ETH 3 

ETH4 

GR 

MH 

MHE 

MHH 

MHZ 

NA 


Exponential integral constants 

Absorption coefficient, cm -1 

Exponential integral constants 

Mass fraction of H 
Mass fraction of H 

e 

Mass fraction of H£ 

Mass fraction of H£+ 

Boltzmann constant, erg/°K 

DENSITY, g/cm 3 

Free-stream density, g/cm 3 

Planck's constant, erg-sec 

Exponential integrals of order 1, 2, and 3 

Defined in Appendix Al 

Euler's constant 
Mole fraction of H 
Mole fraction of H 

e 

Mole fraction of H, 

Mole fraction of H-,+ 

Number density of H ? , raolecule/cm ! 


NA1 

P 

PH 

PHT 

PREI 

QO 

QR 

ROH 

S 

T 

TS 

U 

UI 

V 

VELO 

VF 

VI 

Y 

YD 

YI 

Z 

ZETA 


Number density of H 2 +, oolecule/cm 3 
Pressure, dyne/cm 2 
Sensible enthalpy, erg/g 
Total sensible enthalpy, erg/g 
Free-stream pressure 

Radiative heat flux at shock front, erg/cm 3 -sec 
Radiative heat flux at any field point, erg/cm 2 -sec 
Absorption cross section, cm 2 

Nondimensional coordinate measured along the shock surface 
Temperature, °K 
Shock temperature, 9 K 

Velocity component tangent to the shock surface, cm/sec 

Free-stream velocity component tangent to the shock surface, 
cm/ sec 

Velocity component normal to the shock surface, cm/sec 
Free-stream velocity, cm/sec 
Frequency, sec -1 

Free-stream velocity component normal to the shock surface, 
cm/ sec 

Coordinate measured normal to the shock surface, cm 

Photodissociation yield 

Photoionization yield 

Altitude, km 

Shock angle 
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APPENDIX B1 


PROGRAM "PERC" - TO COMPUTE PRECURSOR EFFECTS 
USING PERTURBATION METHOD 


PROGRAM PERC. (INPUT* OUTPUT * T -'Pc5 = I NPUT * TAPE6 =OUTPUT ) 

DIMENSION 3ATA(4)« 0(4). AA (4) « AR(4). AC(4)« AD ( 4 ) ♦ P2 ( 4 ♦ 2 1 « P2 ( A » 

12)* V ( 4 « 2 ) « VI (4*2)* PC ! O ( 4 ) • D C I 1 ( 4 ) * GO (4)* 0 1 ( A 1 * T ^(4)i T 1 ( 4 ) « 

0T(4)* E ( 8 0 ) * Y I ( 4 1 « Y O ( 4 ) « POE (4 !i Ai(4)« 81(A). Cl (4li r l Ul< C? ( 

1 4 J < o 1 < 4 ) • 03 ( A ) 4 0 0 ( 4 ) < A F < 4 ) « A - ( A ) • A G ( 4 ) • A H ( 4 ) « A ^ ( 4 ) « A I O ( 4 1 « E X ( 7 « O 1 • 

1 SUV ( 6 ) « I ERR (3 ) .01 1 (4 ) ♦ Z 1 ( lO ) 

REAL VHH . VH2 * VHE ♦ '-'H 
EXTERNAL FX1 * FX2 
zt ( 1 1 = 0 .001 
Z1 (2 ) =C * 002 
Z1 (3 ) =C .C9r> 

Z 1 ( 4 ) = 0 • 0 1 

Z1 (5 1=0.02 
Z1 <6 1=0.05 
Z1 (7 1 =C • 1 
Z1 (8 1 =C .2 
Z1 (O ) =0 .5 
Z 1 (10 1 = 1.0 
U = 8 • 3 1 4 3E 7 
TH=6.6256F-27 
GARMr 1 .433 
ROE (1 1 =4 . 1 E- 1 8 
RC^ (2 1 = 6 • 2 E - 1 8 
RCE (3 )=2* 1 E-16 
DC 4 K = 1 .27 

READ (5.15) DEN2. V I NF . T3« PR I « QO 

1 5 FORMAT ( 5E 1 1 • 2 1 

GN^2 = T. 24 3 I ! 2* PR I /! AF • 1 0.** 16*0 . 89 

T ( I ) =RCE ( 1 1 ♦Gnm2 
7 ' ) = ROE ( 2 1 *3NH2 

T " 1 -ROE ( 3 1 *GNH2 
RC AP=25. 

A0=0 .26777343 
A 1 =Q.63A 7608925 
A2= 1 H. 0590 1 6973- 
A3=6. 5733287401 
R 0 = 3 . 956469228 
51 -2 1 .0996530... 

82=25 .632956 1 4 
83=9.5733223454 
DO 8F 7 1=1.1 0 
Z = Z1 (VI >*2 K .* (-1 . 1 
77 IF ( AtiS ( Z 1 .LE * 1 .0F-« ) GO TO 
DO 11 N = 1 . .3 


3-1 


; l 3N*l l ^3f(Wrv-y-w 


«rc$fcWj I, 


<-■ 

l- 

*f : 


I 

? 


i 


I 


F 


| 


i 


X=-Z*T < N ) 

IF (X.LT. 1 . ) GO TO 13 


E < 1 ) = EXP ( -X ) * 
1 *x**3+X**4 ) ) 


EX (N, 1 )=E ( 1 ) 

GO TO 14 

13 E ( 1 > =-0 .5772 1 5664-ALOG (X)+X-X**2./4.+X**3./l 8 • -X**4 • /96 • + X** c . /6O0 
1 . -X**6. /432C .+X**7./35230.-X**8 • / 322 560 . + X**9 • / 326 592? • -X** 1 O ./ 362 
288000. 


ex <N. i > =e ( : > 

14 DO 28 1=1,6 

E ( 1 + 1 ) = (EXP ( -X )-X*E ( I ) ) /I 
EX (N« l + l )=E ( 1 + 1 ) 

28 CONTINUE 
1 1 CONT I NUE 
GO TO 501 
555 EX< 1 , 2 ) =1 . 

EX <2,2 ) =1 . 

EX ( 3 * 2 ) = 1 . 

EX ( l . 3 ) =0 .5 
EX (2,3 ) =0.5 
EX ( 3 , 3 ) =0.5 

50 1 GMARC = V INF/ (70.*(145.**0.5)*100. ) 

P=1 . 385E- 1 6 

V ( 1 , 1 ) =5.02E+1 5 
V ( 1 ,2 ) =8 . 70E+ 1 5 
V ( 2 , 1 ) = 3 • 75E + 1 5 
V (2 , ? ) =5.02F+ 1 5 
V ( 3 , ! ) = 1 • 15E+15 
V ( 3 , 2 ) = 3 • 75E+ 1 5 

VI <1 ,1 >=V(1*1 > *6.6256F-27/ (I .3805^-16*75) 
VI (1 .2 ) =V ( 1 . 2 ) *6 • 6256E —27/ ( 1 .3605F-1 6*TS ) 
VI (2,1 ) = V ( 2. 1 ) *6.62562-27/ (1 .38n?(=_j6*TS) 
VI ( 2 » 2 ) ='/ < 2 * 2 ) *6 • 6256E-27/ (1 .3d05E-16*TS> 
VI (3,1 ) = V < 3 » 1 ) *6.62562-27/ ( 1 .3806E- 1 6*TS ) 
VI <3»2)=V<3,2) *6.62562-2 7/ ( 1 .3805E- 1 6*TS ) 
XI =0 , 

X2=C. 

X 3 = 0 . 

X4 = 0. 

X5 = 0 
X6 = ' n 

X 7r f) 

X8 = 0 


I 
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X9=0 
X IC=C 

DO 500 N= l . 3 
EPS =0.001 

CALL ROMSS (VI CN.l ) « V l ( N » 2 ) « FX 1 » E P S ♦ SUM ( N ) • I ERR ( N ) ) 

CALL ROMBS (VI (N. l ; .vl (N .2 ) ,FX2»EPS.SUf-MN + 3 ) , IERR (N+3 > ) 

'a/RITE (6.32) SlJM(N) .Sl.'M (\' + 3 ) 

32 FORMAT ( 3X . *P2»* .5 1 0. 3. 3X . *P?=*« El o.3) 

TO = 1 4 5. 

GHC= 1 . 527*'J*T0 
Wl =3 .230E-24 
G 1 =1 4.868EI 2 
GD=4 .5E 1 2/2 . 

V I ( 1 ) = 1 . C 
Y I (2 ) =0 .875 
Y I ( 3 ) =0 .0 
YD ( 1 ) =0 . 

YD (2 )=C . 1 25 
YD< 3 ) = 1 .0 

A I (N ) =GNH2*RCE (N)/(DEN2*VINF*GiiO ) 

B I <N ) =- (C,NH?*R0E (M ) * VI )/ (PEN2*V ! NF*GH0 ) * ( ( G I -0 .89* U* TO ) * Y 1 ( N ) + ( 2 
1GD-1 » 8 9 *U * TO ) * Y D ( N ) )/2. 

02 = 1 5/ ( 2. 1 A 1 59**4 )*00 

D1 (N ) =EX ( N. 3 ) *S'JY (\ ) 

02 (N ) =YD ( N ) *FX <N, 3 ) *8t.iM ( N+3 ) / ( P* IS ' 

C3 <N ) =Y I ( N ) *EX (N. 8 )*SIJM ( M + 3 ) / ( P* TS ) 

D 1 1 (N>=Q2*D(N) 

X 1 =D (N ) +X1 
X 2 = D 1 ( N ) +X2 
X 3 = X3+D2(N) 

X4 =X4+D3 ( N ) 

500 CONTINUE 

BA T = 2 • / (0.7?7*GARM*GVARC**2-1 . ) 

VZ=-02*BAT*X 1 
PRF=2. *GAPM*C2*X 1 
DEN=-VZ 

GH T = 2 • / < DEN2 * V I NF ** 3 ) *02*X2 
CH2= 1 • *'//! *2. *02* X4/ ( v I NF*OFM2 ) 

C H = £ • * '*' 1 * 0 2 * X 3 / ( V ! NF *DEN2 ) 

R = 0 . 

VP = r> • 

8 1 ' TF 1 = 1 45 •* (PRE-DEN- (Ch2*O)*20C!./'180.I7> 

CHH = .: .801R-O-2-CH 



C HE = 0 • 1981 

r e l 1 1 1 =TF1 +14?). 

RT=Tp1 111*1 .98726 

IF ( TF I 1 I 1 .GT. 1 000. ) GO TO 999 

IF (TF 1 1 1 1 .GT.30C. ) GO TO 998 

EThP? = CH2 *( ( TF 1 1 1 1-1 OC • ) /20G • * ( 1 23? • 1+358305. 1/2. 

ETHP) =CHH* ( ( TF i 1 1 1 - 1 nc . )/?90.* < 1 26/1 .81-1 264 .81/2. 

E T HP 3= CH*( (TF1 11 1-100. )/2C!0 . * ( 966 . 1 + C 1 1 34 . ) 

ET+P4 = CHE* ( ( TF1 1 1 1-1 00. )/2CC.* (984.4 1-984 .41/4. 

GO TO 1001 

998 ETHP! =CHH*RT* (2.8460849+4. 19321 1 6E-3*TFl 111/2.-9.611 9332E-6+TFJ 1 1 1 

1 **2/3»-9.51 52974F— 9*TF 1 1 1 1 **3/4 • +3 • 3 09 3492F- 1 2*Tp 1 1 1 1 **4/=; .-Q. 6 7^6 
2372E2/TF1 1 1 1 1/2. 

ETHP2=CH2*RT* (2.81 7375+3. 6"761E-7+TF1 ’ 1 1 /2 . -7 . 0*6* 4PF-6*T F I 1 1 1**7/ 
1 3.+9.26140E-9*TFl 1 1 1 **3/4 .-3 .09C228E- 1 2* TF 1 1 1 1 **4/S.+l .80C273E6/TF 

2 11111/2. 

GO TO 1000 

999 ETHP 1 = C.HH*RT* (3.0436897+6. 1 1871 1 0E-4*TF1 1 1 1 /2 . -7 . 3993 c 5 1 F -9* T- 1 1 1 1 
1 **2/3.-2.0 33 1 90 7~- 1 1 *TF 1111 **3/« . +2 . 45 a 3 70 j =•- \ 5*TP 1 1 1 1 **4/5 .-8.^49 
21 0-2E2/TF1 111 )/2. 

£T-P? = CH2*RT* ( 3.328 71 5+2. 605067E-3* TP 1 1 1 1 /2. + 1 • 42 2 452 OF — 7* TF l 1 1 1 ** 
12/3.-4.459C24E-1 1 *TF 1 1 1 1**3/ 4. +3. 7 33^56E- 1 6*TF 1 1 1 l**a/ c .+l .7997470 
2E5/TFJ 111 )/2 . 

1 uOG ETHP3 = CH*RT* (2.E + 2.5473497E4/TF1 1 1 1 1 

PTHP4 =CH5*RT* (2.6-7. 4 5374 F2/TF1 1 1 1 )/4. 

1 00 1 £TH = ETHP1 +ETHP2+PTHP.3+ETHD4 
ETH<=ETH*4 .184*1000.0 
WRITE < 6 . 161VINF. DEN2.TS 

16 FORMAT (//4X.*VINF=*,E1 1 .3. 4X. *OEN=* . E 11 . 3 . 3 X . *T S=* . E 1 1 .4 ) 

■VR I TE (6 *52 1 

52 FORMAT < 8 X.*R*. 1 OX .*Z*. 8 X.*VZ*. 8 X.+PRE*. 8X.+0EN*. 8 X.+T-+. Bx , *0. 

1H2 + *. 7 X . *CH * . 7X.*GHT*. 8 X.*ETH+) 

itfR I TE ( 6 » 37 )R . Z* VZ .PRE . DEN * TF1 ♦ CH2 . CH ♦ GHT , E TH< , Z 1 (Ml ) 

37 FORMAT (/2X.2F1 0.2.9E1 1 .3 ) 

GO TO 80 

Z = 7+ 0,10 

IF (7.GT.+.91 j GO TO 4 
GO TO 77 
68 CONTINUE 
4 CONTINUE 

38 STOP 
ENP 

F'JNCT ION c X 1 ( v ) 

PX 1 •= /** 3/ ( ExP ( X 1 -1 1 


return 

FNO 

FUNCT I ON FX2 ( x } 
p;<? = X * +2/ (EXP ( X 1 - 1 1 

R p T i |R N 
FNO 


l 


t 

I: 
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APPENDIX B2 





PROGRAM "THIN" - TO COMPUTE PRECURSOR EFFECTS 
USING THIN-LAYER APPROXIMATION 


PROGRAM Thin ( input, OUTPUT. TAPE5= INPUT, TAPE6=0UTduT ) 

COMMON ARF<3 ) .SUM ( 12 ) . I 1 « T ( 50 ) ♦ TS . V I ( 3 ) . Y0 < 3 ) • NA ♦ ROH ( 3 I ♦ I ERR < 1 2 ) 
I ♦ VF<3»2 ) ♦ OH * CK iC'<' OC « Y , S » Z3 ( 3 ) .JiVELO I , DFN I 
DIMENSION V (50 ) »U(5C ) «P (50 ) .PHT ( 50 > »PH(50 > .CH2 (50 ) »CH(50 ) .DFN(5r> ) 
REAL NA.Nr5.NAl .mh2.vhh.mh.vhE 
CO 2 1 N= 1 ,20 

READ (5,1) Y.S.TS.PREI , Z T A , Z » 00 , VELO I .DENI 
FORMAT (4F8.2.2F5.3.E1O.3.E0.O.Z1O.3 ) 

VELO I =-VELO I 

ZTA = ZTA*3.1 415026535/(2. *90. ) 

ROH ( ! ) =4 • 1 E-l 9 
ROH ( 2 ) =9.25-1 8 
ROH ( 3 ) =2 • 1 E- 1 8 
DH=6.6256E-27 
C<= 1 • 38C54E- 1 6 
C=3.CE 1 0 
VF (1,1) =8 . 7E 1 5 
VF(2,2)=3.75E15 
VF ( 3 , 1 ) =3.75515 
VF (3,2 ) =1 . 15E1 5 
VP( I .2 ) =5 . 02E 1 5 
VF ( 2 , 1 ) =5 . 02E 1 5 

V I = VEL 0 I * S I N ( Z T A ) 

U I =Vr LO I *C05 ( ZT A • 

V < 1 ) = V I 
U ( 1 ) = U I 

P ( I ) =PRE I 
DEN ( 1 >=DENI 
CH? (l ) =0. 

CH ( l ) =0 . 

A = 1 4 . 8E 1 2 
0 = 4. 5E 1 2/2 . 
vi ( 1 )= 1 .0 

V I (2 ) =0 .875 

VI ( 3 ) =0 . 

YQ( 1 ) = C .0 

YO ( ? ) =0 • 1 25 
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Y0(3 ) = 1 . 

R=8.34 1 3E7 

PHTtl ) = ( VELO I ** 2 ) /2 • + 1 .5?7*R*1 45 • 

Pw( 1 )=1 .537*144. *R 
T ( 1 ) = 1 45. 

20 DO 1 1 1 1=1*49 

11 = 1 

NA = ( 7.243 1 1 22E22*C.89 )*P(1 )/(T(l )*10.) 

NA 1 = NA* 1 . OE-6 
ARF ( 1 ) = RCH ( 1 ) *NA 1 
ARF(2 )=RQH(2 1 *NA 1 
ARE ( 3 ) =ROH ( 3 ) *NA 1 
DEN ( 1+1 >=DENl*Vl/V( I ) 

V ( I + 1 > = ( DEN I * ( V 1 **2 ) -P ( l ) +P < 1 ) ) / < DEN I *V 1 ) 

U ( I + 1 > =U I 

CALL QRAC1 A < Y,S«0P) 

dht I s 1 .52 _, *p* 1 45. + ( VELO I **2 )/2* 

PHT ( 1 + 1 )= (DENI *VI *PHT! -OR)/ (DENI *VI ) 

OH ( I +1 ) sPHT ( I + 1 ) - (V ( I + l )*+2+U(I+l )**2)/2. 

T( i + j ) s ( oh ( 1 + 1 1— (5./4.*R+T( I )+A/2. )*CH2( I )-(3./4.*R + T( I ) +D ) * CH ( I ) ) 
1/(1 .527 *R ) 
call PCH2 (PCHI .PCHD ) 

NB = 1 . 

CH ( I + 1 ) =N5*PCH0 
CH2 ( I + 1 ) =N8* = CH I 

□ (j + 1 ) = DEN ( I + 1 ) *R*T (! + !)*( 18-0. 17/4 CO. +o.4*(CH2(I + l ) +CH ( I ♦ 1 )) ) 

Jr (ABS( ( V t 1+1 ) — V < I ) ) / V ( I ) ).(jl .'.J.Ul ) ou 

IF ( AE S ( ( T ( I + 1 > — T ( I ) 1 / 7 ( I ) ) • GT • 0 • 0 l ) GO TO 59 

IF ( ASS ( ( P ( I + 1 ) -R ( I ) ) /P ( I ) ) . 3 T .0 • •'l ) GO TO 59 

IF ( AaS ( (PH ( I + 1 ) -oh ( I) )/ch ( I )) .GT ,C .0 1 > GO TO 40 

IE ( AaS ( (PHT ( I -t-1 ) -PHT ( I ) ) /PhT ( I ) ) .GT .C .C 1 ) GO To ^9 

IE { ABS < (CH( 1 + 1 J-CH ( I ) )/CH ( I ) ) .GT.O* 31 > GO TO 59 

I f ( AES < < CH2 ( I + 1 ) -CH2 < I ) 1 /CH2 ( I ) ) .GT «C .0 1 ) GO TO =9 

I F ( AbS ( < DEN ( I + 1 ) -DEN ! I ) ) /DEN ( I ) ) .0 T . G . 0 } ) GC TO 5 9 

GO TO 1C 
=9 GO TO 111 
1 1 1 CONT I N'JE 
1C WRITE (6*11) S.Y 

11 format (4X. +S=* *E1 0 .3 » x **CY* .5/ *+Y=* ,E 1 0 .3 . X . +cv* ) 

WRITE (6*12) Z. VELO!. CD 

1 2 FORMA T ( c X . * A TT I T ODE = * « •- 1 C • 3 » 4 X . * VEL CO I T Y I \|F = * . E 1 D . 3 « * C v / SC * . 4 X . + 

1 9 ( 0 5 =* . El C. 2 * * ER-G /C v 3 — LED * . ✓ / ) 

CHH = ... • MO 1 R-CH2 ( I *■ 1 ) -CH { I + ’ ) 

CHE=0 .1961 


V 
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AD0»CHH/2.+CH2 ( l + l )/2.+CH( I+l )+CHE/4. 

MHH* ( CHH/2 • )/AOD 
MH2«(CH2( l+l )/2« ) /ADD 
MH«(CH( l+l ) )/ADD 
MHE» (CHE/4 l/ADO 
RT*1 .98726*T( l+l ) 

IF (T ( 1+1 ) «GT. 1 000. ) GO TO 999 
IF < T C l+l ) .GT *300 • ) GO TO 998 
TFlll 1»T( 1+1 ) 

ETHP2»MH2*< (TFl 1 1 1 - 1 00 • ) /20C.* ( l 232 • )+3583P5. ) 

ETHPl =MHH*( (TFlll l-l 00. >/200.*( 1264*8 )-1264.8 ) 

ETHP3»MH*< < TFl l t 1-100. )/200.*(966. >+51 134, ) 

ETHP4 = MHE*< (TFl 1 l 1 - 1 OC • ) /200.* (984 .4 ) -984 .4 ) 

GO TO l 00 l 

998 ETHPl *MHH*RT* (2.8460849+4. 1932,1 16E-3*T( I +1 ) /2 • -9.6 1 1 9332E-6*T (l+l ) 

1 **2/3* “9 .51 22974E-9*T ( 1+1 )**3/4.+3.3C9C492E-l 2*T ( 1+1 > **4/5. -9.6725 
2372E2/T ( 1+1 ) ) / 

ETHP2»MH2#RT*(2.817375+3.65761E-3*T( 1+1 )/2.-7.965546E-6*T ( l+l )**2/ 
1 3.+S.26 1 40E-9*T ( 1+1 > **3/4 .-3.09C228E-1 2*T ( f+l )**4/5.+l .800273E5/1"( 
21 + 1 ) > 

GO TO 1000 

999 ETHPl *MHH#RT* (3.0436897+6. 1 ie7110E-4*T( 1 + 1 > /2 . -7. 399355 l E-9*T ( l + 1 ) 
i**2/3.-2.0331 907E-1 1 *T ( I + 1 ) **3/4 .+2 .459379 1 E“1 5*T ( 1+1 )**4/5.-8.549 
2 1 002E2/T ( I+l ) ) 

ETHP2»MH2*RT* (3.328715+2.505067E-3*T( I+l )/2. + l .4224520E-7*T (l+l >** 
12/3.-4.459024E-1 1 *T ( I + l )**3/4.+3.' , 33756E-15*T ( I + l )**4/5.+l .7997470 
2E5/T ( l+l) ) 

1000 ETHP3®MH*RT* ( 2 • 5 t- 2 • 54 70497E4/T ( I + ! ) ) 

ETHP4*MHE*RT*(2.5-7.4 53T4P2/T( I + l ) ) 

1 001 ETHs (ETHPl +ETHP2+ETHP3+ETHP4 )*4. ! B4* 1 0000 ./2 . 3 1 4 
WRITE (6 »28 ) ETH 

28 FORMAT (4X»*ENTH0PY ( ABSQL ) = * » El 1 .4*/) 

WRITE (6,13) 

1 3 FORMAT (3X«*V*» 1 l X»*U*» 1 2X » *T* , 10X.*P* , 1 C X . *DEN* • °X » *H* . 1 CX»*HT* , l 1 
1X.*CH*. ICX,*CH2*) 

WRITE (6»14) V ( I + l ) ,U( I + l ) ,T( I +1 ) .P ( I + l ) , DEM ( I+l ) , PH ( J + j ) .PHT ( I *1 ) 
1 , CH ( 1 + 1 ) «CH2 ( I + l ) »QR , I 

14 FORMAT (X,El0.4,2X,E9.3,2X.2X,E9.3,2X,E9.3,2X,EIO.3,2x,El'',3,2X.El 
l0.3,2X,E10.3,2X,El0.3,3X.E11.3,3X.F3.t,///) 

IF (Y.GT.0.5) GO TO 15 
V*Y+0. 1 
GO TO 20 

15 Y® V+C • 4 

IF (Y.GT.2. ) GO TO 21 
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GO TO 20 
21 CONTINUE 
STOP 
END 

SUBROUTINE OR AO 1 A ( Y *S »QR ) 

COMMON ARF ( 3 ) » SUM ( l 2 ) . I . T ( 50 ) . TS • V I (3 > « yq(3 > .NA.R0H(3> •IE»R(12)«VF 
1 (3*2) » OH « C< *C«<I »QO»YY«SS»tJ (3) 

COMMON /FFF/ Z 

DIMENSION E2 <3 ) ♦ ARY (3 > *E 1 ( 3 ) . VI ( 3 %? ) 

EXTERNAL FX 1 « FX2 • FX3 
Fl*6.256E-27/( 1 •3805E-16*TS ) 

VI (1*1 ) = VF ( 1 .1 ) *F 1 
VI (1 .2 ) = VF ( 1 « 2 ) *F 1 
VI (2.1 ) =VF (2*1 ) *F 1 
VI ( 2 • 2 > = VF (2.2) *F 1 
VI (3.1 ) = VF ( 3 « l ) *F 1 
VI (3*2) * VF (3.2) #F 1 • 

EPSa0.002 
DO 50 <=1.3 
<1 =K 

CALL ROMBS ( VF (< ♦ 2 ) « VF (<»l),* r Xl« E°S • SUM ( < ) * I ERR ( < ) ) 

IF ( IERR<< > .EO.O) GO TO 30 
PRINT 31 .<. I ERR ( < ) 

30 CALL ROMBS ( VI (< • 2 ) . VI (< . 1 ) »FX2 . fcPS « SUM (<+3 ) . 1 fcWK ( K+3 ) ) 

IF ( IERRK + 3 ).EQ.O) GO TO 40 

<3=<+3 

PRINT 31. <3. lE»R(<+3) 

40 A=0.0001 
Z=V+0*0001 
e=v 

IF (B.EG.O. ) GO TO 32 

CALL ROMBS (A.B. FX3.EPS. SUM <k+6> . I ERR ( <+6 ) ) 

IF ( 1ERR(K+6).EG.J) GO TC 50 
< 6=<+6 

PRINT 31. <6. IER«(<+6 > 

31 FORMAT (*<=*. 12 .* IERR=»« 13 ) 

GO TO 50 

32 SUM ( <+6 ) * 0 • 

50 CONTINUE 

OR 1*0 

00 I 00 < = 1.3 

GR=0.5772 

ARY(<)=A«F(<)*Y 

IF ( Y.LT • C • CO 1 ) GC to 33 


GO TO 20 
21 CONTINUE 
STOP 
END 

SUBROUTINE OR AO 1 A ( Y *S »QR ) 

COMMON ARF ( 3 ) » SUM ( l 2 ) . I . T ( 50 ) . TS • V I (3 > « yq(3 > .NA.R0H(3> •IE»R(12)«VF 
1 (3*2) » OH « C< *C«<I »QO»YY«SS»tJ (3) 

COMMON /FFF/ Z 

DIMENSION E2 <3 ) ♦ ARY (3 > *E 1 ( 3 ) . VI ( 3 %? ) 

EXTERNAL FX 1 « FX2 • FX3 
Fl*6.256E-27/( 1 •3805E-16*TS ) 

VI (1*1 ) = VF ( 1 .1 ) *F 1 
VI (1 .2 ) = VF ( 1 « 2 ) *F 1 
VI (2.1 ) =VF (2*1 ) *F 1 
VI ( 2 • 2 > = VF (2.2) *F 1 
VI (3.1 ) = VF ( 3 « l ) *F 1 
VI (3*2) * VF (3.2) #F 1 • 

EPSa0.002 
DO 50 <=1.3 
<1 =K 

CALL ROMBS ( VF (< ♦ 2 ) « VF (<»l),* r Xl« E°S • SUM ( < ) * I ERR ( < ) ) 

IF ( IERR<< > .EO.O) GO TO 30 
PRINT 31 .<. I ERR ( < ) 

30 CALL ROMBS ( VI (< • 2 ) . VI (< . 1 ) »FX2 . fcPS « SUM (<+3 ) . 1 fcWK ( K+3 ) ) 

IF ( IERRK + 3 ).EQ.O) GO TO 40 

<3=<+3 

PRINT 31. <3. lE»R(<+3) 

40 A=0.0001 
Z=V+0*0001 
e=v 

IF (B.EG.O. ) GO TO 32 

CALL ROMBS (A.B. FX3.EPS. SUM <k+6> . I ERR ( <+6 ) ) 

IF ( 1ERR(K+6).EG.J) GO TC 50 
< 6=<+6 

PRINT 31. <6. IER«(<+6 > 

31 FORMAT (*<=*. 12 .* IERR=»« 13 ) 

GO TO 50 

32 SUM ( <+6 ) * 0 • 

50 CONTINUE 

OR 1*0 

00 I 00 < = 1.3 

GR=0.5772 

ARY(<)=A«F(<)*Y 

IF ( Y.LT • C • CO 1 ) GC to 33 


AO-0.26777343 

Al *8.6347608925 

A2*l 8*0590169730 

A3*8. 5733287401 

80*3 *958469228 

81*21.09965308 

82*25*6329561 A 

83*9.5733223454 

IF C ARY (K)*UT*1*) GO TO 200 

X*APY(K ) 

El CK )*EXP(-X )# CAO+A 1 *X+A2*X#*2+A3*X**3-fX**4 ) / <X* t tJ U+tS l *X+fc52*X**2+d 
13*X**3+X**4 ) ) 

GO TO 220 
200 X*ARY ( K ) 

El CK ) *-0*57721 566 -AlOG ( X )+X-x**4 */*4 • t***j./ 

1 •— X**6*/4320 * + X**7*/3528C • -X**3 */322560 * + X#*9 */3265920* — X**l 0./362 
288000* 

220 E2CK >* CEXPC-X )-X*El (K) ) 

E3 CK )* CtXPC— X )— X>F2 CK ) )/?• 

GO TO 35 
33 E2CK )*1 • 

E3CK) = 0.5 

35 GR2*2** ARFC<)*DH/<C**2* >*SUMCK+6 ) 

QR4*E3 C K ) #SuM (K+3 ) * 1 5 • *00/ (3*141 59*»'o ) 

ORl *OR2+OR4+OR 1 
100 CONTINUE 

QRsQRl *2. *3* 14159 

RETURN 

ENO 

FUNCTION FXt (X ) 

COMMON ARFC3 ) « SUM ( 1 2 ) « I » T C 50 ) » T 5 t Y I ( 3 ) » YO ( 3 ) »NA*UOH (3 ) ♦ IFR* t 1 2 > « VF 
1 C3»2)«OH»CK*C.K«00 
FX1 =X**3*EXP C- <DH*X/ <C«#T ( I) ) ) ) 

RETURN 

ENO 

FUNCTION FX2CX) 

COMMON AKFC3 ) «SUM C 12 ) « l «T(50)»TS.YI <3 > • YDC 3 > «NA ♦ RC* <3 > « I CkK C 1 2 ) » VF 
1 C 3 « 2 ) • OH • C< ♦ C « K » Q 0 
0AT*5*6697E-5 
FX2»X**3/CEXPCX )-l • 1 
RETURN 
END 

FUNCTION FX3CX) 

COMMON ARF ( 2 1 • SUM ( 1 2 1 « I • T C 50 ) • TS • Y I (3 > • YD C 3 > *NA,*.:C H C 3 ) . ! CPU ( 1 ? ) « VF 
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1 ( 3 « 2 ) • Dh . CK » C . < . Q 0 
COMMON /FFF / Z 

FX3- ( 1 •+ ( • 577-1 • + ALOG< ARF ( < )*<Z-X ) ) )* ( ARF <<) * (Z-X ) j — ( ARF (< )* (Z-X 
1 ) ) *#2Z2*+ < ARF (< ) * < Z-X ))**3/ 12* ) *SUM (< ) 

RETURN 

END 

SUBROUTINE PCH2 (PCHI.PCHD) 

COMMON ARF < 3 ) .SUM ( 1 Z )♦ J » I ( 5U )♦ * P» * 1 » J )« YU i J ) «N« »HOin j )• 1 LH* i 1 2 ) »V* 
1 (3*2) .DH.C<»C«K.Q3.Y.S.E3(3 ) »J1 .VELOI .CENI 
D I MENS I ON VI (3.2) .0-11 (3) »C.H2 (3» 
external FX4 

FI =6.256E-27,’ ( 1 .3805F-16*TS ) 

VI (1 . 1 >=VF< 1 ♦ 1 )*F1 
VI U .2)=VF( 1 .2>*F1 
VI (2.1 ) = VF ( 2 ♦ 1 )*F1 
VI (2.2) =VF < 2 .2 ) *F 1 
VI (3.1 )=VF(3.1 ) *F 1 
VI (3.2 ) = VF (3.2 ) * c 1 
EPS=0«0C1 
DO 59 J=1 .3 
J1 =J 

CALL ROMBS ( VI ( J.2 ) . VI ( J« 1 ) ♦ FX4. EPS. SUM ( J+9 > ♦ .* ERR ( J+9 ) ) 

CHI ( J ) sY I (U ) *E3 ( J )*SUM ( J + 9 )/ ( l ,3QO=5E-l 6*TS ) 

CH2( J) *YD ( J )*E3 ( J >*Sl iM ( J+o ) / ( 1 * 33n5E — 1 6*TS ) 

59 CONTINUE 

PDHi=CHl ( 1 ) +CH 1 (2 )+CHl (3 ) 

PDHD*CH2 ( 1 ) +CH2 (2 )+CH2 ( 3 ) 

X 1 *GG* 1 5* *2. *3 .23E-24/ ( DEN I * VELD 1*3. 14159**4 )* (- 1 . ) 

PCH I sPQH | *x 1 
PCHOsPDHO*Xl 

return 

END 

FUNCTION FX4 ( X ) 

COMMON ARE ( 3 ) . SUM (12 ) . 1 . T ( 50 ) * T S . Y 1 (3 J . YD ( 3 J «NA . ROH ( 3 ) ♦ I - RR ( 1 2 ) . V p 
1 ( 3 » 2 ) » DH . C< ♦ C ♦ < • 0 3 • Y . 5 . E3 ( 3 ) .J 
FX4=X**2/(EXP(X)-1 ) 

RE TURN 
END 
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